Large-scale dynamical simulations of the three-dimensional XY spin glass 
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Large-scale simulations have been performed in the current-driven three-dimensional XY spin 
glass with resistively-shunted junction dynamics for sample sizes up to 64 3 . It is observed that 
the linear resistivity at low temperatures tends to zero, providing a strong evidence of a finite 
temperature phase-coherence (i.e. spin-glass) transition. Dynamical scaling analysis demonstrates 
that a perfect collapse of current-voltage data can be achieved. The obtained critical exponents agree 
with those in equilibrium Monte Carlo simulations, and are compatible with those observed in various 
experiments on high-T c cuprate superconductors. It is suggested that the spin and the chirality order 
simultaneously. A genuine continuous depinning transition is found at zero temperature. For low 
temperature creep motion, critical exponents are evaluated, and a non-Arrhenius creep motion is 
observed in the low temperature ordered phase. It is proposed that the XY spin glass gives an 
effective description of the transport properties in high-T c superconductors with d-wave symmetry. 

PACS numbers: 05.10.-a, 68.35.Rh, 74.25.-q 



I. INTRODUCTION 

The ordering of the three-dimensional (3D) XY spin 
glass model [lj has attracted considerable attentions. 
Theoretically, earlier work suggest transitions at zero or 
very low temperatures. Following the pioneering work by 
Villain 2] where the role of chirality (i.e. vorticity) arising 
for non-collinear XY system, was emphasized, Kawamura 
et al proposed a chiral-glass transition, even though the 
finite temperature spin-glass transition vanishes Q . Later 
on, by calculating the Binder ratio and the spin-overlap 
distribution function for lattice sizes up to L — 16, Kawa- 
mura and Li have given a numerical evidence that the low 
temperature phase is a chiral glass without the conven- 
tional spin glass order Q. 

However, this spin-chirality decoupling scenario is 
not consistent with several studies @, B 0, B IS 
Maucourt and Grempel[5| and subsequently, Akino and 
KosterlitzQ found evidence for finite temperature spin- 
glass transition from zero-temperature domain wall cal- 
culations. By performing a finite-size scaling analysis of 
the correlation length £, the most successful technique 
developed in the Ising spin glass 11 1 , Lee and Young Q 
observed a transition at the same temperature for both 
spins and chiralities for lattice sizes up to L = 12. Beside 
the equilibrium simulation, the resistive behavior in the 
XY spin glass has been also studied. In a vortex represen- 
tation, Wengle and Young showed evidence of a resistive 
transition at finite temperatures [12J for lattice sizes up to 
L = 10. Granato observed that the current-voltage char- 
acteristic in the XY spin glass with both bimodal and 
Gaussian couplings for the lattice size L — 12 exhibited 
scaling behaviour^, [13] , which was interpreted in terms 
of both the phase-coherence (spin-glass) transition and 
the chiral-glass transition. 

The controversy on the spin-chirality decoupling sce- 
nario is however still ongoing [r|. For the Heisenberg 



spin glass, the transitions at the same temperature was 
also observed for both spins and chiralities in Ref. Q 
with the use of the same method for XY spin glass. 
However, subsequently, it is found that the situation 
becomes rather unclear when results of larger sizes are 
included 0, EUGl G3, EI- The data at the low temper- 
atures in larger samples show rather marginal behavior, 
i.e. the system is close to the lower critical dimension. 
Recently, there seems a intensely competition for the lat- 
tice sizes accessible 0, [lj], Ell, the record until now 
is L = 48 18]. Motivated from the debate for Heisen- 
gurg case, Young and his collaborator further study the 
XY spin glass with larger lattice sizes up to L = 24(l9j. 
and observed similar marginal behavior for low temper- 
atures and large systems. So it is also desirable to per- 
form the simulations on very large sample for XY spin 
glass. In contrast with the Ising and Heisenberg spin 
glass, there is an additional technique in XY spin glass, 
that is the measurement of the resistivity from dynamical 
simulations [tj E3, El]- Moreover, one can access larger 
system in this kind of dynamical simulations, because 
a steady state is more easily reached than equilibration 
duo to the presence of the external driven force. It is also 
of interest to know whether the previous picture based 
on small systems {L < 12) in the dynamical simulation 
, El] is modified on larger systems for the XY case. 

On the one hand, the XY spin glass has found ex- 
perimental realization, not only in layered manganite, 
e.g. EuQ.^Sr^MnOi [2(J, but also in high-T c cuprate 
superconductors with the d-wave paring symmetry [2l| . 
which can be regarded a random distribution of 7r 
junctions [22l [23l . uA , H^, [26j|. The nature of d-wave sym- 
metry will changes the sign of the coupling between XY 
spins, while the spin angle denotes the phase of the su- 
perconducting order parameters. So the XY spin glass 
is expected to be used to interpret the phenomena ob- 
served in high-T c cuprate superconductors, such as vor- 



tex glass phase [13, |H, [2^|, a true superconducting state 
with vanishing linear resistivity. The evidences to sup- 
port the existence of this phase have been reported in 
many experiments (30l. l3l| . 

Since the chiral variable in the XY spin glass can be 
defined as the vortex in the plaauette[9l[To|. some tech- 
niques developed in superconducting vortex models (3^. 
|33| . |34| can be in turn employed to characterize the na- 
ture of the low temperature ordered state in the XY spin 
glass. As is well known, in the random pinning envi- 
ronment, the energy landscape for the vortex motion is 
highly nontrivial. The theoretical understanding for the 
nonlinear dynamics response has been advanced in many 
years [35j . But the full theoretical study is still very chal- 
lenging, computer simulations are hopeful to provide use- 
ful insights. 

In this paper, by resistively-shunted-junction dynam- 
ics, we perform large-scale dynamical simulations in the 
3D spin glass, both the phase-coherence transition tem- 
perature T g and the critical exponents are estimated. 
The depinning transition at zero-temperature and creep 
motion below T g are also investigated. The rest of the 
paper is organized as follows. Sec. II describes the model 
and dynamic method. Sec. Ill presents our main results, 
where some discussions are also performed. Finally, a 
short summary is given in the last section. 




FIG. 1: Current-resistivity curves at various temperatures 



In the xy plane, the supercurrent between sites i and j 
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(A Xl A y ) the fluctuating twist variable and 
Oi = fa + rj • A. The new phase angle 9{ is periodic in 
both x- and y-directions. Dynamics of A„ can be then 
written as 



A « = h E K 



+ -la^OL = x,y (3) 
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II. MODEL AND DYNAMIC METHOD 

The Hamiltonian of the 3D XY spin glass in the phase 
representation is given by Q 



H 



cos 



{fa -fa)- 



(1) 



where the sum is over all nearest neighbor pairs on a 3D 
square lattice, fa specifies the phase of the superconduct- 
ing order parameter on grain i, denotes the strength of 
Josephson coupling between neighboring grains with zero 
mean and standard deviation unit. The present simula- 
tions are performed with the lattice size L = 64 for all 
directions, considerably larger than those in literature. 

The Resistivity-Shunted- Junction dynamics is incorpo- 
rated in simulations, which can be described as 



ah 
~2e 



dH 



Je: 



E 



Vij > 



(2) 



The voltage drop is V = —LA X . 

The above equations can be solved efficiently by a 
pseudo-spectral algorithm [32j due to the periodicity of 
phase in all directions. The time stepping is done using a 
second-order Runge-Kutta scheme with At = 0.05. Our 
runs are typically (4 — 8) x 10 7 time steps and the latter 
half time steps are for the measurements. Our results are 
based on one realization of disorder. The present system 
size is much too larger than those reported in literature, 
a good self-averaging effect is expected. We have per- 
formed an additional simulations with a different realiza- 
tion of disorder for further confirmations, and observed 
quantitatively the same behavior. Note that it is practi- 
cally hard to perform any serious disorder averaging for 
the present very large system size. 



III. SIMULATION RESULTS AND 
DISCUSSIONS 

A. Finite temperature phase-coherence transition 



where J ex t,j is the external current which vanishes ex- 
cept for the boundary sites. The rjij is the thermal noise 
current with zero mean and a correlator (t)Vij (t')} = 
2akBTS(t — t'). In the following, the units are taken of 
2e = h = a = ks = 1. 

In the present simulations, a uniform external current 
I x along rr-direction is fed into the system, the fluctuating 
twist boundary condition [34j | is applied in the xy plane, 
and the periodic boundary condition is used in the z axis. 



First, we study the phase-coherence transition. The 
current-voltage characteristics are simulated at various 
temperatures ranged from [0.2, 0.6]. At each tempera- 
ture, we try to probe the system at a current as low as 
possible. The voltage is determined when a steady state 
is arrived. Fig. 1 demonstrates the resistivity R = V/I 
as a function of current / at various temperatures. It 
is clear that, at lower temperatures, R tends to zero 
as the current decreases, which follows that there is a 
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FIG. 2: Dynamic scaling of current- voltage data at various 
temperatures according to Eq. (UJ). 



FIG. 3: Dynamic scaling of current- voltage data at T g — 0.2 
according to Eq. ©. 



true superconducting phase with zero linear resistivity. 
While R tends to a finite value at higher temperatures, 
corresponding to an Ohmic resistivity. These observa- 
tions provide a strong evidence of occurrence of a phase- 
coherence transition at finite temperature in the 3D XY 
spin glass in the dynamical sense. 

For a continuous phase transition characterized by the 
divergence of the characteristic length and time scales 
t ~ £ z (z is the dynamic exponent) , Fisher, Fisher, and 
Huse [28[ proposed the following dynamic scaling ansatz 

TRZ z+2 - d = tf ± (Jf , - 1 /T). (4) 

where d is the dimension of the system (d = 3 in this pa- 
per), and £ oc| T/T g — 1 \~ u is the correlation length which 
diverges at the transition, ^(x) is a scaling function, with 
+ and - signs corresponding to T > T g and T < T g . Eq. 
(QJ was often used to scale measured current-voltage data 
experimentally [30l |3~H . 

To extract the critical behavior from the numerical re- 
sults of the current-voltage characteristics, we will also 
perform a dynamical scaling analysis. As shown in Fig. 
2, using T g = 0.33±0.02, z = 4.0±0.1, and v = 1.4±0.1, 
an excellent collapse is achieved according to Eq. Q. 
The value of v is by no means close to the Ising value 
2.15±0.15[ll|, demonstrating the XY and Ising spin glass 
belong to different universality classes. 

The finite-size effect is particularly significant at tem- 
peratures sufficiently close to T g when the correlation 
length exceeds the system size. For the temperatures 
considered and the very lattice size L = 64 here, we be- 
lieve that the finite-size effect is negligible in the present 
simulations. To confirm this point, we perform partic- 
ular simulations right at T g = 0.33 obtained above for 
different system size. At T g , the correlation length is cut 
off by the system size in any finite system, the scaling 
form ((3]) for d — 3 becomes 

TgRL*- 1 = y(IL 2 /T g ). (5) 

As shown in Fig. 3, a good collapse is shown using 
z = 4.2 ± 0.1. This consistence demonstrates that the 



estimate from Fig. 2 is reliable. Therefore a new evi- 
dence of a finite-temperature phase-coherence transition 
is provided convincingly in the 3D XY spin glass. 

The above obtained T g and static exponent v are con- 
sistent with T g = 0.34 ± 0.02 and v = 1.2 ± 0.2 obtained 
in equilibrium Monte Carlo simulations |8[ for lattice sizes 
up to L = 12. More recently, the data for the spin-glass 
and chiral-glass correlation length at larger sizes show 
marginal behavior, indicating that the lower critical di- 
mension is close to 3 for both spins and chiralities[l9j. 
The crossing temperatures for the spin and chiral-glass 
correlation length decrease very slightly with the size 
with large uncertainties, close to 0.3 for largest lattice 
sizes they can simulate. The data for the ratio of the 
chiral and spin-glass correlation for the largest sizes in- 
tersects for T about 0.33 and then splay out in the low 
temperature side. The size dependence of the crossing 
temperature at large sizes is very weak, supporting a sin- 
gle transition for both spins and chiralities. Interestingly, 
this crossing temperature is in excellent agreement with 
the present T g for a very large lattice size. We do not 
think it is a coincidence. 

Previous driven Monte Carlo dynamical simulation on 
the same model for a small lattice size (L = 12) [9( esti- 
mated T g — 0.335 , v = 1.2, and z — 4.4, basically consis- 
tent with the present ones. It follows that the results are 
not sensitive to the system size in the dynamical simula- 
tions, even the detailed dynamics is different. The Monte 
Carlo dynamical simulations in the vortex representation 
also show an equilibrium transitive transition 12]. The 
estimated value for v — 1.3 ± 0.3 agrees with the present 
one. But the dynamic exponent z « 3.1 is obvious lower, 
indicating different dynamic universality class for the XY 
spin glass in the phase and vortex representations. 

Typical values of the correlation length exponent and 
the dynamical exponent extracted from the various ex- 
periments on high-T c cuprate superconductors [30l l3~l| 
fall in the range of v « 1.0 — 1.7 and z w 4.0 — 5.0. 
Our results for these static and dynamic exponents in 
the XY spin glass are consistent with these experimen- 
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FIG. 4: Log-log plots of V — I curve at zero-temperatures. 

tal ones. In addition, nonlinear resistivity measure- 
ment in the Y Ba2Cu±0% bulk sample [24j| near the on- 
set of the paramagnetic Meissner effect have been in- 
terpreted as a chiral-glass transition owing to the pres- 
ence of 7r junctions. The nonlinear contribution p<i to 
the resistivity was observed to have a peak at the tran- 
sition with power-law behavior p ~ J~ a and exponent 
a — 1.1 ± 0.6, which has been reproduced by a gen- 
eralized XY spin glass model[26j]. Granato proposed 
that it could also be interpreted as a consequence of 
the underlining phase coherence transition and the expo- 
nent a is determined by the dynamic exponent through 
a = (5 — z)/20. By the above obtained dynamic expo- 
nent, we have a = 0.5 ±0.05, compatible with the exper- 
imental value within the error bar. Note that the present 
value of z in the large-scale simulations is smaller than 
those previously given by Granato[l, [T^ by about 10%, 
so our estimated a is more close to experimental one. 
Based on the above comparisons with experiments, we 
propose that the XY spin glass gives an effective descrip- 
tion of the transport properties near the phase-coherence 
transition in d-wave superconductors. 

Until now, at least there is a general consensus in 3D 
XY spin glass that the chiral-glass transition occurs at 
a finite temperature larger than or equal to the spin- 
glass transition temperature T g 0, i, B, 0, !, H, E3, til- 
The data for the chiral correlation length in the re- 
cent large-scale equilibrium Monte Carlo simulations 
reveals that the chiral-glass transition temperature is 
around 0.30 — 0.33, not larger than the present spin-glass 
transition temperature. We therefore suggest that the 
spin and the chirality probably order simultaneously. 



B. Depinning transition and creep motion 

With the low temperature spin-glass phase in hand, 
in the remaining part of the paper, we will study the 
depinning and creep phenomena in this phase. 

To study the depinning transition at zero tempera- 
ture, we start from high currents with random initial 
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FIG. 5: Log- log plots of V — T curves at three currents around 

phase configurations. The current is then lowered step 
by step. Fig. 2] shows the current- voltage characteristics 
at T = 0. Interestingly, we observe continuous depinning 
transitions with unique depinning currents, which can be 
described as V oc (7 - I^f with I c = 0.115 ± 0.001, 
(3 = 1.99 ± 0.02. Note that the depinning exponent (3 is 
greater than 1, consistent with the mean field studies on 
charge density wave models (36j. 

At low temperatures, the current-voltage characteris- 
tics is rounded near the zero-temperature critical current 
due to thermal fluctuations. An obvious crossover be- 
tween the depinning and creep motion can be observed 
around I c at the lowest accessible temperatures. In order 
to address the thermal rounding of the depinning tran- 
sition, Fisher [3^ first suggested to map this system to 
the ferromagnet in fields where the second-order phase 
transitions occur. This mapping was latter extended 
to the random- field Ising model [33], flux lines in type-II 
superconductors [Hj], and Josephson junction arrays [38). 
If the voltage is identified as the order parameter, the 
current and temperature are identified as the inverse tem- 
perature and the field in the ferromagnetic system respec- 
tively, analogous to the second-order phase transitions, a 
scaling relation among the voltage, current, and temper- 
ature reads 

V(T, I) = T^^iT- 1 ^ 5 ^ - I c /I)l (6) 

where S(x) is a scaling function. 

Eq. (6) reveals that right at I c the voltage shows 
a power-law behavior V(T,I = I c ) cx T 1 / 5 , providing 
a method to determine the critical exponent 1/6. The 
V — T curves at three currents are presented in Fig. [5] 
on a log-log plot. The critical current is seen to be be- 
tween 0.11 and 0.125. The values of voltage at other cur- 
rents within (0.11,0.125) can be evaluated by quadratic 
interpolation. The square deviations from the power 
law can be calculated. The critical current is defined 
at which the square deviation is minimum. We obtain 
I c = 0.115 ± 0.02, consistent with those obtained above 
at zero temperature. The temperature dependence of 
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FIG. 6: Scaling plot of the current-voltage data at various 
temperatures below T g according to Eq. ©. 



voltage at the critical current is also plotted in Fig. 
yielding 1/8 = 1.315 ± 0.001. 

With the critical exponent 8 and the critical current 
I c , we can adjust the depinning exponent [3 to achieve 
the best data collapse according to the scaling relation 
Eq. (JBJ) for I < I c . In Fig. [6l a optimal data collapse 
of the current-voltage data at various temperatures be- 
low T g provides an estimate of/3 = 2.00±0.01, which is 
in excellent agreement with those derived at T = de- 
pinning transition. Moreover, the scaling function with 
the form V oc T 1 / 6 exp[A(l - ^/T 05 } is used to fit well 
the current-voltage data in the creep regime, which is 
also list in the legends of Fig. [6l As we know, the 
product of the two exponents /35 describes the tempera- 
ture dependence of the creeping law. 08 « 1.52 deviates 
from unity, demonstrating that the creep law is a non- 
Arrhenius type. The non-Arrhenius type creep behav- 
iors have been also previously observed in charge density 
waves [39j. 

In the 3D XY spin glass, the existence of a stable 
spin glass phase at finite temperature is well established 
through both previous equilibrium studies [H, [l|| and 
present dynamical simulations. The value of (36 is close 
to 3/2, similar to the 3D vortex glass and gauge glass 
model[34[. We believe that the value of (38 is generally 
3/2 in the glass phases in various strongly disordered 
3D frustrated XY model. The further analytical work is 
highly called for. 



IV. CONCLUSIONS 

We have performed large scale dynamical simulations 
on the 3D XY spin glass within the resistively-shuntcd- 



junction dynamics. The strong evidence for the low tem- 
perature spin-glass phase is provided in dynamical sense. 
By the dynamical scaling analysis, a perfect collapse 
of simulated current-voltage data is achieved by using 
T g = 0.33 ± 0.02 , z = 4.0 ± 0.1, and v = 1.4 ± 0.1. 
These critical values agree with those in the previous 
equilibrium Monte Carlo simulations as well as driven 
Monte Carlo dynamical simulations. In the resistive sim- 
ulations, the size dependence of various exponents and 
critical temperature in the phase-coherence transition are 
very weak at larger size, providing a good technique to 
measure critical properties in this model. The static and 
dynamic exponents are compatible with those observed in 
the various experiments on high-T c cuprate superconduc- 
tors, which underlines the significance of XY spin glass 
in d-wave superconductors. Combined with the recent 
equilibrium Monte Carlo simulations, we suggest a single 
glass transition involving both spins and chiralitics. 

We also study the depinning transition at zero temper- 
ature and creep motion at low temperatures in detail. A 
genuine continuous depinning transition is observed and 
the depinning exponent is evaluated. With the notion of 
scaling, the critical exponents are estimated, which are 
consistent with those from independent simulations at 
zero temperature and at critical current. The value of 
(36 is close to 3/2 and the scaling curve is fitted well by 
a exponential function, suggesting a non-Arrhenius type 
creep motion in the low temperature ordered phase in 
the XY spin glass. 

Finally, it is proposed that the XY spin glass model 
may capture the essential transport feature in high-T c 
cuprate superconductors with d-wave symmetry Further 
experimental and theoretical works are clearly motivated. 
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